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Abstract 

We analyze a one-component simple fluid in a liquid- vapor coexistence 
state, which forms an arbitrarily curved interface. By using an approach 
based on density functional theory, we obtain an exact and simple expres- 
sion for the grand potential at the level of mean field approximation that 
depends on the density profile and the short-range interaction potential. 
By introducing the step-function approximation for the density profile, 
and using general geometric arguments, we expand the grand potential in 
powers of the principal curvatures of the surface and find consistency with 
the Helfrich phenomenological model in the second order approximation. 
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1 Introduction 



Although the description of a one-component fluid in a liquid- vapor coexistence 
state has been studied since long ago, it still remains as a topic of interest for 
inhomogeneous fluids [2 [H [3J |3] . Not just open questions are yet to be answered 
within this problem, but also its appropriate description may be useful for the 
analysis of more complex systems such as colloids, polymers, mixtures of these, 
etc. Particularly, investigation of the microscopic expression for the free energy 
that represents an arbitrarily curved interface remains a topic of interest for 
the one-component fluid in a coexistence state. In order to analyze this sys- 
tem different approximation schemes, rigorous derivations, and many genuine 
shortcuts have been implemented [11 [5l El [71 [8] . Associated to the free energy 
there is the calculation of microscopic expressions for interfacial properties such 
as surface tension, spontaneous curvature, and rigidity coefficients. The main 
objective of this work is to derive the free energy and all these quantities in a 
rigorous and fully general form. 

One of the approximation schemes that has provided a first principles descrip- 
tion is the stress tensor theory, implemented initially by Romero- Rochin, Varea, 
and Robledo ^9^. In this theory the authors identify the normal component of 
the stress tensor as the fundamental quantity to calculate the grand potential, 
which represents the energetic cost to maintain the interface having a given ge- 
ometry. Construction of the most general form of this stress tensor is due to 
Percus and Romero-Rochfn, and represents one of the important achievements 
within this approximation scheme |10[ . The microscopic expression for this 
stress tensor has been used to describe the interfacial region, within the van 
der Waals model, when the interface has planar, spherical, and cylindrical ge- 
ometries. For each of these interfaces, the exact microscopic expression for the 
grand potential has been constructed considering drops of arbitrary size |12] . 
Nevertheless, in order to compare these results, it is convenient to approximate 
the grand potential as a limiting case; which is, when the radius of the Gibbs 
dividing surface is much larger than the range of the interaction potential. This 
has lead to a microscopic expression for the grand potential as a power series of 
the inverse radii of curvature. At this level, we establish an equivalence between 
the microscopic expression for the grand potential and the Helfrich phenomeno- 
logical model 13 for each of the corresponding geometries. 

Despite the Helfrich model is widely used to describe a variety of systems [3 
I14j . no first-principles derivation of the model is known. Several efforts have 
been made aimed at achieving the task [T5]. The usual strategy to get to the 
phenomenological form of the Helfrich free energy is based on the assumption 
that the interface is a bidimensional elastic continuous-medium [16] . According 
to this model, the free energy of the interface is a function of the principal 
curvatures of the system, which can be written as 




(1) 
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where H = (l/Ri) + (I/-R2) and K — I/R1R2 are the mean and Gaussian cur- 
vature respectively, withi?i and R2 being the principal radii of curvature. The 
coefficients 7, cq, k, and R are equilibrium surface properties. 

As a second part of this work, we carry out a rigorous derivation of the Hel- 
frich model corresponding to a fluid interface, Eq. ([T]), calculate the surface 
properties, and show that it is a legitimate model to describe arbitrarily curved 
interfaces. We follow an essentially different program as compared to previous 
approaches. Instead of analyzing diverse particular geometries, which is not 
a simple task, we consider the general case of an arbitrarily curved interfacial 
region and study it using a semi-orthogonal coordinate system. A notable fact 
that appears in the description of curved interfaces concerns the localization of 
the Gibbs dividing surface. As the physical properties of the system must be 
independent of this choice, a displacement in the localization of the interface 
in microscopic distances shall not modify the value of the surface tension, but 
will alter the value of the rigidity coefficients. Thus, different localizations of 
the Gibbs dividing surface give rise to different values of the rigidity coefficients, 
which results in an arbitrariness of their values. Explanation of the origin of this 
arbitrariness is one of the questions yet to be addressed. In this study, we fix the 
radius of the Gibbs dividing surface within the approximation introduced for 
the density profile. Although this criterion is not unique, the results obtained 
are consistent with previous works. 

This paper is organized as follows. In Section [2] we briefly outline the stress 
tensor formalism that describes a smooth, but otherwise general, interface. Sec- 
tion [3] is devoted to the calculation of the microscopic expression for the grand 
potential of a curved interface within the van der Waals approximation. Next, 
the limit of large radii of curvature, as compared to the range of the interac- 
tion potential, is introduced in Section |4] and a further approximation to the 
grand potential is obtained, which coincides with the Hclfrich phenomenological 
model. Finally, some concluding remarks are drawn in Section [51 



2 Stress Tensor Theory 

According to density functional theory, the grand potential can be written in 
the form [2 [2] 



where F[p{'r)] is the intrinsic Helmholtz free energy, fj, is the chemical potential, 
and Vext is the external potential. The equilibrium value of the density profile 
is obtained by minimizing the grand potential; that is by solving equation 



n[p{r)] = F[p{r)] + / df\p - Kxt(r)]p(f), 



(2) 




Po 
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In absence of exact analytic solutions, the usual approach is through numerical 
methods, which provide estimate values. However, in this work we insist in an 
analytic solution and follow the route of the stress tensor to obtain it. Next we 
briefly outline the general formalism, which can be read in detail in Ref. [5]. 

In an inhomogeneous fluid the condition for mechanical equilibrium implies the 
existence of a conservation equation, obtained from ([3]) , which may be expressed 
as 

V-a = poVycxt, (4) 
with a being the stress tensor and jOoK:xt(^) the external force per unit area. 

The stress tensor is not unique. As it can be seen from Q, a term with zero 
divergence may always be added. Nature of the system suggests a separation of 
the stress tensor in two contributions 



CT = 0-0 + 0-5, (5) 

where uo is the homogeneous part and as is the inhomogeneous one. The ho- 
mogeneous part describes the bulk phases, where the density profile is uniform, 
whereas the inhomogeneous region is that where \/po{r) 7^ 0. 



The free energy of the system is obtained by integrating the normal component 
of the stress tensor over the whole space. Once again, we separate this normal 
component in two contributions, one associated to the homogeneous region and 
the other from the inhomogeneous one [9] 

^[Po{r)] = / dra^'ir) = / dfa^ (r) - f dra^{r). (6) 



As we are interested in the inhomogeneous region, we neglect the contribution 
to the energy arising from the bulk phases. From now on we concentrate in 
obtaining the microscopic expression for fig, which is given by 

^^5bo(r)] = - /dra^(r-). (7) 



We observe that is the fundamental quantity needed to obtain microscopic 
expressions for the properties of the surface. For the system under study the 
microscopic expression for the stress tensor of the inhomogeneous region, within 
the van der Waals approximation, is [TTJ [H] 

= - J dr' j\xpo{r-{l~X)r')Coi\f'\y^Vppo{r + Xr') 

- iv.y dr' j\\po{T-{l-\)r')C:{\T'\) 

X r'p[ry,po{r + \r')-ryc,po{r + \r')]. (8) 
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This quantity depends exclusively on the density profile, the interaction poten- 
tial, and is independent of the geometry of the interface. In fact, the geometry is 
defined by the density profile. For instance, for a planar interface po{r) = po{z), 
and for a spherical interface po{r) = po(l^l)- In general, for an arbitrary interface 
if ^ denotes the normal coordinate, the density profile is a function exclusively of 
this quantity: po{r) = Po{0- By using this information, we go on to calculating 
the grand potential for an interface having an arbitrarily curved geometry, 



3 Free Energy of a Curved Interface 

The free energy of the intcrfacial region, given by Eq. ([7]), in explicit form writes 



(9) 



- ^ j drh^{f)fip{r)V, J dr' J dXpo{r - {1 - X)r')Q{\r'\) 
X r'i^KV,po{f+ Ar') - ry^poif+ A?-')]. 



By performing some manipulations this expression may be compacted to the 
form [H] 

rJs = i I dr I dr' I dXp(){f-{l-X)r')i:j{\f'\)VaPa{r + Xr')[r'i^na{r)np{r) 



2Xr'pnp{r)r'^V afiuir) + Xr'pna{r)r'^V ^fip{r) + Xr'pnfi{r)r'^V ^ha{r)], 



(10) 



which has been previously used to describe interfaces having planar, spherical, 
and cylindrical geometry |12| . Direct calculation of the grand potential in each 
of these geometries leads to the equation 



= I dr I df' 



ds Vpoif) ■ Vpoif')wis + {f- f'f), (11) 



the difference being contained only in the density profile. 



The first goal in this paper is to prove that Eq. ((TT|) is satisfied for an arbi- 
trarily curved interface, and we concentrate on this for the remaining of this 
section. Instead of starting from Eq. (jlOp . we prefer to manipulate expression 
([9|) , introduce the linear change of variables 



'+Ar', 



f(2) ^ ^' 



and use the relationship 



j,(2)v(l) 



Af(2) 



d_ 
dX 



(12) 
(13) 



(14) 
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with the consideration of the symmetry with respect to the exchange of indices 
a and (3 in this equation, and that the profile is constant at ±00, to obtain 

Q.S = ~ j df j df'po{j^w{\f'\)ha{f)n0{f)r'0'^aPo{r + f') 



dr J d.rj dXXw{\r'\)r'i^r'^ha{r)np{f) 
X V^po(r - (1 - X)f')\7^po{r + Ar ') 
- ^/'^^/^'"/ dXXw{\-r'\)r'i^r'^ha{f)nis{f) 

X po(r-(l-A)f=")V.V.po(r + Ar="). (15) 

The last two integrals in this expression have to be manipulated to eliminate 
the arbitrary parameter A. With this in mind we define a system of semi- 
orthogonal coordinates, having unit vectors {h,ti,t2), and express each of the 
vectors f and f' on its own basis. That is, r = rnn{r) + rt-^ti + ri.^t2 and 
r ' = rn'n'{r') + rfj^x + rt^i2- However, vector r ' can also be expressed on the 
basis of vector r in the following manner: r ' = (r ' • h{r))h(r) + (r ' • ti)ti + (r ' • 
^2)^2 = '''nfi'i^ + ^ti'^i + ''/,2^2- We now go on to manipulating expressions so as 
to eliminate A. To achieve it we define an auxiliary function W{r'^ + r'^^ + ^'ta); 
related to il;(|r^^ + rf^ + ^'^a I) via 

irW^T'^ + r'l +r[l) = r'^w{\r'^ + r',\ + \). (16) 

n 

This relationship is employed to eliminate any power of within the integrals. 
After some more manipulations we get to the different expression for the grand 
potential 



= dr j dr' po{f)w{\r'\)ha{r)hi3{r)r'pVaPo{r + r') 

- \j dfj dr' y ' dAA dt w{t + r'^ + r',\ + r',l) 
X WuPo{r-{l-X)r')S/^po{r + Xr') 

- ^JdrJdr'J dXxJ dtw{t + r'^ + r't^) 
X po(r- (1 - A)r')V^V^po(r + Ar') 

- ^J^^J<^^'J '^^^ j j dt'w{t + t' + r'^ +r'tl + 
X ^[V-/'o(r - (1 - A)f ')V,po(r + Af ')] 



i J dr J dr' J^' dXX d.t dt' w{t + t/ + r'^ + r'^) 

[po(r - (1 - A)r")V.V.po(r + Ar')], (17) 
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dr!? 
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which may now be simplified by considering each of the terms separately. As 
it is shown in Appendix \^ the last four integrals cancel one another in pairs, 
yielding the result for the grand potential 

^S = -^ J drj df'poir)w{\r'\)ha,ir )npiryp\7^po{f+ f '). (18) 

The task now is to show the equivalence between this and Eq. (fTT|) . We start 
by carrying out further manipulations using (jl6p and introducing a change of 
variables of the form (fT2)) -([T3 | . but now being 

f(i) = f+r', (19) 
f(^) = f. (20) 

Thus, the final result obtained for the microscopic expression for the grand 
potential of an arbitrarily curved interface is 

ns = -\j dr^^^ J df^'^ p dtw{t + (r^i) - ri2))2 + _ fi^)f) 

X 9iiVo(r(i))ai^Vo(r1^)), (21) 

where each vector has been written as r ''^ = rn^v}-^^ + rt^if' + rt^i!^^ ■ 

Notice that this expression is equal to pip , which was to be proved. This is one 
of the most general results for the grand potential that represents the free energy 
of the interfacial region. The result is simple and exact. Given the interaction 
potential, it only depends on the density profile. 



4 Local Approximation on the Free Energy 

Although the previous result is simple and exact, it is not the appropriate form 
to carry out comparisons with other works. In this sense, it results convenient 
to get an approximation as a reference value. We start by choosing a common 
basis to express all vectors. If the basis is that of vector r'^^\ vector r^^^ 
can be written as f (2) = (f (2) . ^(i))^(i) + (^(2) . + (f (2) . ^(i))^^ ) = 

ria^n^^^ + r-t^i^^T^ + ^-^^^^T^- By introducing these into ^ we find 

1 nOO p pOO p poo 

ns = ^ dri^^ / dS^''> / dr(2) / / dt 

X w{t + (r,W - r . ^(i))2 (^-a) _ ^^(2))2) 
X nW.n(2)9(i)po(f(i))arWpo(r(''-n(i)), (22) 

where drl^2 — d/drn] and each element of volume has been written as dr^"^^ = 
dvn^dS^^^ with r/*'' two-dimensional vectors on 5'-^' and S^'^^ respectively. 
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The expression for the grand potential that we have constructed so far contains 
the description of the whole interfacial region, independently of its details. The 
microscopic expression is exact for this model; exact in the sense that the den- 
sity profile contains all the information of the geometry under consideration. In 
order to obtain microscopic expressions for interfacial properties it is necessary 
to fix the Gibbs dividing surface. The choice in this work is through the approx- 
imation of the density profile as a step function. For a more general profile the 
procedure becomes more complicated. We will consider that case in a future 
study. 

The density profile is thus approximated as 

Po(r) - Po(r;«) = po;e(r(:) - r^) + po.e(r« - r«), (23) 
where r^g is the radius of the Gibbs dividing surface. We have 

5«Po(r«)--Apo'5(r«-r«), (24) 

with Apo — pQy — Pol. Evaluation of the derivatives of the density profile is 
obvious. The result we obtain is 

ns = \{Apofj dS^'^ I d5(2)n(i).n(^)^°°di^(i+(r«-r(fj)^+(rl«-r^^^^^ 

(25) 

To proceed we approximate the surface representing the interface as a paraboloid. 
This is possible as long as the radius of this surface is very large compared to 
the range of the interaction potential. Under this assumption, the surface can 
be approximated as a plane with corrections |17[ 118) . 



We now build a local coordinate system in the neighborhood of a point P of the 
Gibbs dividing surface. This arbitrary point is localized on the surface by vector 
f^^\ The normal vector at that point is n*-^-* and we chose the z axis pointing 
in this direction. Thus the coordinates x and y lie on the plane tangent to the 
surface at P and point along the directions of the principal radii of curvature. 
At point P, f/^^ = (0,0) and n'^-* = k, so that the Gibbs dividing surface is 
localized at height r^Q*. Vector r'^^ localizes another point Q on the dividing 
surface, close to P but outside the local plane, at a distance z — rj^Q — rljo* from 
the plane (see Fig. [T]). If z is measured from the local system, it has the value 



2 Y i?i i?2 

As vector r is not parallel to the z axis of the coordinate system, the normal 
vector at point Q is given by 



( 2 V- 1 

V -Ri ' -R2 ' ^ 



(27) 
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Figure 1: This schematic picture shows the local approximation for the surface 
S about a point P. Points P and Q are localized by vectors f and r '^^^ 
respectively, and the normal vectors to the surface at those points are n^^^ and 
h^'^\ Point P is chosen as the origin of the local coordinate system whereas Q is 
outside the tangent plane. Its projection onto this tangent plane has coordinates 
{x,y). The radius of localization of the Gibbs dividing surface is at r^o'*. The 
distance from Q to the local plane (projection of n^"^^ onto the direction n*-"^-*) is 
r^o — = z, with z seen from the local system as a paraboloid. 
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The metric in this coordinate system is 5 = 1 + [\7z{x, y)]^ = 1 + ^ + ^ i where 
<C 1 and -C 1. The surface element in the local system is dS^'^'' — g^dxdy. 

The scalar product of the normals is fS^^ ■ h'-'^^ = . We incorporate the effect 
of the local approximation into the interaction potential to obtain the free energy 
of the interfacial region 



By expanding the interaction potential to first order about < + + and 
evaluating the integrals on the coordinates x and y, one gets to a microscopic 
expression for the grand potential in terms of the principal radii of curvature 



00 Jo 
2 „,2 ^ 2 



ns = - dSl'-^^ / drCo{r^)r 



where we have put dS'*^^^ = dS. On the other hand, from the definitions of the 
mean {H) and Gaussian [K) curvatures one finds that 



^H^--.K = is2 + i^ + ^nnr- (30) 



4 1 1 

-K = — H 

3 i?J i?2 3i?ii?2 

By introducing this into ((29)) . we get to the most general microscopic expression 
for the grand potential within this approximation level; consistent with the 
Helfrich prediction [13] for the free energy 



j dS 






[ 2 7o 



64 V 3 

(31) 

which is also in agreement with previous results in the planar [51 [71 [T^j , spheri- 
cal [H], and cylindrical [12] geometries, as it is shown in Appendix iB] 



5 Concluding Remarks 

There exist two relevant aspects of this work that are worthwhile remarking. 
The first one concerns a rigorous proof for a general expression, exact and sim- 
ple, of the grand potential within a mean field approximation. The second is a 
first principles derivation of the Helfrich free energy, within this context, from 
a completely original approach. This, in addition, confirms that the Helfrich 
scheme is appropriate for the study of curved interfaces. Although the local 
approximation of a surface as a plane that has been used is appropriate for the 
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description of weakly curved surfaces, we observe that the result for the grand 
potential representing the free energy of the system is sufficiently general. The 
microscopic expression obtained is a function of the principal curvatures of the 
surface, in complete agreement with previous predictions (13[ [T5] . We also find 
complete consistency with previous results for the simplest geometries, which 
were obtained using a different analytic approach |T2]. Finally, we point out 
that within the step function approximation for the density profile, no contri- 
bution exists from spontaneous curvature to the free energy. We shall study the 
problem of an arbitrarily curved interface for a more general density profile in 
a future publication. 
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A Simplifications on the Grand Potential 

In this appendix we summarize the main steps to simplify Eq. (|17p . We start 
by writing it in the form 

ns = J df J dr'po{r)w{\r'\)na{r)ni3{f)r'ijVaPo{r + r 

where the T,;, with i = 1, . . . , 4, have been defined as the integrals appearing in 
PT|) sequentially. 

In order to simplify and easily identify each of the terms Ti, T2, T3, and T4, we 
introduce explicitly the change of variables (|T2t - (fT3|) . Under such a coordinate 
transformation the term Ti reads 

Ti^-\j df^'^ J dr^'^ y'"dtc:;(t+|f(2)|2)a(i)po(f(i)_f(2))9(i)p„(f(i))^ 

(33) 

where the notation dji^ — d/drn^ has been introduced. Integrating by parts 
with respect to rn \ this now becomes 

Ti=i ydf(i) ydf(2) j\xx y"dic^(<-Kir(2)|2)p„(f(i)_f(2))a(i)2^g(^-(i)). 

(34) 
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The other term with normal derivatives of the same order is T2. The effect of 
transformation (fT2l)-(IT3l) into this leads to 



J Jo Jo 

(35) 



4 , 



Therefore these terms cancel one another when substituted into the expression 
for the grand potential ([5^ . 

To simplify the terms with normal derivatives of the second order, that is T3 
and T4, we need to take into account the following manipulations 



d< d< dr^l^ d< Sri') 
from where 

92 d 



= A- 



d 



dr 



(1) 



+ -^ = A9,W+9(2), (36) 



dr[^ dr. 
By introducing this into T^, we obtain 
1 



^ (Aa(i) + a(2)) = A2a(i)2 + 2A9,wa(2) + d^^)\ (37) 



= -1 I dr^'> I dr^'> I dXX I dt I dt'£l;(t + t' + |f(2)|2 



00 poo 

dt 



A2a(i)2 + 2Aa(i)a(2) + ai^)^] b,Wpo(r-(i) - f (2))aa)po(r-(i))l (38) 



We need to calculate derivatives of the product within the integrand. The first 
one is 



+ 

and the second 



5a)2^^(,^(i)_,-(2))a(i)^^(,-(i)) 
aWpo(rW-r(2))aWV„(,-(i))(39) 



Cna(^''po(r(^)-r(2))9(i)p„(f(i)) 
9,Wpo(r(i)-f(2))9(i)2po(f(i))] 
9,W3po(r(i)-r'(2))a(i)po(f(i)) 
2a(i)V„(f(i)_f(2))a(i)2p„(,-(i)) 

CVo(r(i)-f(2))9(i)3p„(f(i)).(40) 



The corresponding derivatives for dn'' are also calculated. By introducing all 
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these into T3 one finds 

T3 = dr^^n dr^^n dXX dt dt'cD(t + t'+|f(2)|2) 

+ a(iVo(r - r (i))] + 2A[a(i)Vo(r - f (2))a,Wpo(f (^)) 

(41) 

or alternatively, after an integration by parts respect to in each term, 

-| /• p pi pOO pOO 

T3 = ^ dr^^n dr^^n dXX dt (ii'(I>(i + i' + |f(2)|2) 

X {A^ [5(1)^0 (f(^) - r (2))a(i)2po(^-(i))2a(i)po(r - r (2))a(i)3po(r (D) 
+ Po(r - f (2))a(i)Vo(r (^)] + 2A5^)[a^)po(r - f(^))a^)Vo(r (^)) 

Now we investigate the effect of the same transformation in the expression for 
T4. Direct substitution implies 

■\ p p pi pOO pOO 

T4 = ~ dr^^n dr^^U dXX dt cii'a;(i + t' + |f(2)|2) 

X [a25(i)2 + 2AS(i)9(2) ^ 5(2)2j J^^(^-.(i) _ .(2))^a)2p„(f (1))] . (43) 

Although this compact form appears simple, we need to carry out further sim- 
plifications to identify similar terms. We start by calculating the first derivative 

+ po(r(^^-r(2))aWVo(r(^)). (44) 

and then a further derivative of this quantity in the same direction dn'' , which 
yields 

+ Po(r(i)-f(2))aW3po(fW)] 

= a(i)2po(r(i)-f(2))9a)2^^(,-<(i)) 

+ 29(i)po(r(i)-r(2))a(i)3p„(f(i)) 

+ po(r(i)-r'(2))a(i)Vo(r(i)). (45) 
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By substituting we find 



dXX 



dt 



X {A2[9(i)Vo(r - r (2))a(i)Vo(r (1)) + 2di'^po{r^'^ ~ r (2))5(i)3p„(f (D) 
+ Po(r - r (^))5ii)Vo(r (^))] + 2 A^i^) [9(^)^0 (r^^^ - r-(^))9(i)Vo(r (^)) 

= -T3. (46) 
That is, also the terms Ta and T4 cancel one another. 

B Simplest Geometries 

To make contact with previous results we consider here the simplest geometries. 
For the planar surface i?i = i?2 = 00 so that 4iJ^ — |if = and thus Eq. (|3ip 
simplifies to 

'--^[drr^^^r^). 



dS 



(47) 



That is, the only contribution to the free energy comes from the surface tension 
term, in agreement with previous results [5J [7J [T^] . 

For a spherical interface Ri — R2 = R, H ^ j^, and = so that 4iJ^ — 
~ From (PT|) one gets to the corresponding free energy 



1 dS 






[ 2 J, 



Finally, for the cylindrical interface Ri — 00 and R2 — R, so that 4H^ 
In this case the free energy is 



1 dS 






[ 2 io 



drr^ w(r^) 



37r(Apo) 
64i?2 



(48) 



(49) 



These last two expressions are consistent with previous results [12] . 
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